Cosmological data and indications for new physics 
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Data from the Atacama Cosmology Telescope (ACT) and the South Pole Telescope (SPT), com- 
bined with the nine- year data release from the WMAP satellite, provide very precise measurements 
of the cosmic microwave background (CMB) angular anisotropics down to very small angular scales. 
Augmented with measurements from Baryonic Acoustic Oscillations surveys and determinations 
of the Hubble constant, we investigate whether there are indications for new physics beyond a 
Harrison-Zel'dovich model for primordial perturbations and the standard number of relativistic de- 
grees of freedom at primordial recombination. All combinations of datasets point to physics beyond 
the minimal Harrison-Zel'dovich model in the form of either a scalar spectral index different from 
unity or additional relativistic degrees of freedom at recombination (e.g., additional light neutrinos). 
Beyond that, the extended datasets including either ACT or SPT provide very different indications: 
while the extended- ACT (eACT) dataset is perfectly consistent with the predictions of standard 
slow-roll inflation, the extended-SPT (eSPT) dataset prefers a non-power-law scalar spectral index 
with a very large variation with scale of the spectral index. Both eACT and eSPT favor additional 
light degrees of freedom. eACT is consistent with zero neutrino masses, while eSPT favors nonzero 
neutrino masses at more than 95% confidence. 

PACS numbers: 98.80.Es, 98.80.Jk, 95.30.Sf 



I. INTRODUCTION 

A wide variety of observations of Cosmic Microwave 
Background (CMB) and Large Scale Structure (LSS) 
power spectra over the last decade have indicated that 
cosmic structures originated from seed fluctuations in the 
very early universe. The leading theory explaining the 
origin of the cosmological seed perturbations is cosmic in- 
flation [1 j, a period of accelerated expansion at very early 
times. During the inflationary stage, microscopic quan- 
tum fluctuations were stretched to macroscopic scales to 
provide both the initial seeds for the primordial den- 
sity perturbations and tensor (gravitational-wave) fluc- 
tuations [2-6 . Despite the simplicity of the inflationary 
paradigm, the exact mechanism by which cosmological 
perturbations are generated is not yet established. 

In the standard slow-roll inflationary scenario associ- 
ated with the dynamics of a single scalar field (the in- 
flaton), density perturbations are due to fluctuations of 
the inflaton itself as it slowly rolls down along its po- 
tential. In the simplest case, fluctuations are of the adi- 
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abatic type, namely they are sourced by the degree of 
freedom that is dominating the energy density during in- 
flation (the inflaton). In other mechanisms for the gener- 
ation of perturbations, e.g., the curvaton mechanism [7], 
the final adiabatic perturbations are produced from an 
initial isocurvature mode associated with quantum fluc- 
tuations of a light scalar degree of freedom (other than 
the inflaton), whose energy density is negligible during 
inflation. The isocurvature perturbations are then trans- 
formed into adiabatic perturbations when the the extra 
scalar degree of freedom (the curvaton) decays into radi- 
ation after the end of inflation. A precise measurement 
of the spectral index, ns, of the scalar perturbations to- 
gether with a detection of gravity-wave signals in CMB 
anisotropics through its B-mode polarization will provide 
a strong hint in favor of single-field models of inflation. 
Indeed, alternative mechanisms predict an amplitude of 
gravity waves far too small to be detectable by experi- 
ments aimed at observing the £?-mode of the CMB po- 
larization. 

While inflation is the leading candidate model for the 
generation of primordial perturbations, in the words of 
Ref. [8] : "Inflation is at the same time a spectacular phe- 
nomenological success, and an enduring theoretical chal- 
lenge." The phenomenological success is that inflation is 
a simple model for the generation of seed perturbations. 
The theoretical challenge is to understand how inflation 
is embedded in a broader theory or model of fundamental 
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physics. Detailed examination of the CMB perturbations 
are a possible way to discriminate among inflation mod- 
els [9], perhaps even leading to a reconstruction of the 
inflaton potential [TO] . 

Another application of CMB data is to search for ev- 
idence of "new" physics, like additional relativist ic de- 
grees of freedom at recombination [11 , neutrino masses 
[12], early dark energy [13], modified gravity [14], or 
variation of fundamental constants like the fine-structure 
constant [15] or the gravitational constant [T6] . 

The goal of this paper is to examine whether existing 
cosmological datasets can provide evidence for the dy- 
namics of the inflaton field during inflation or evidence 
for new "neutrino" physics. Evidence for the dynamics 
of the inflaton field during inflation would be a depar- 
ture from the Harrison-Zel'dovich (HZ) model (a scalar 
spectral index of unity and no tensor perturbations). 1 
The departure from the HZ model could take the form 
of a scalar spectral index different than unity, a "run- 
ning" (a scale-dependence) of the scalar spectral index, 
or evidence for tensor modes. 

The type of new "neutrino" physics we model would 
be a mass for neutrinos or additional relativistic degrees 
of freedom contributing to the expansion rate around the 
time of recombination. 2 

Our analysis will include CMB data from the nine- 
year data release of the Wilkinson Microwave Anisotropy 
Probe (WMAP9) [19 , the South Pole Telescope (SPT) 
[20] , and the Atacama Cosmology Telescope (ACT) [21], 
including measurements up to a maximum multipole 
number of / max ~ 3000. We will also include informa- 
tion from measurements of baryonic acoustic oscillations 
(BAO) from galaxy surveys in the form of three datasets: 
data release 7 (SDSS-DR7) [22 and data release 9 (SDSS- 
DR9) [23] from the Sloan Digital Sky Survey, and the 
WiggleZ project [24]. We will also use data on the Hub- 
ble constant [25] . 

This study has two motivations. On one hand, the 
Planck collaboration [26] will release soon their first flow 
of data regarding the CMB anisotropics, and therefore it 
is timely to have a state-of-the-art pre-Planck assessment 
of slow-roll inflation. On the other hand, we wish to 
answer three basic questions: 

1. Is the simple Harrison-Zel'dovich model compati- 
ble with current cosmological datasets, or is there 
support for a more complicated perturbation spec- 
trum? 

2. Is standard neutrino physics consistent with cur- 
rent cosmological datasets, or is there support for 



new neutrino physics in the form of a neutrino mass 
in excess of a few tenths of an electron volt or a 
change in the effective number of light neutrinos? 

3. In the event that there is support for physics be- 
yond the HZ spectrum and standard neutrinos, 
can one tell whether the data provides information 
about the primordial perturbation spectrum or the 
neutrino sector? 

As we will see, all combinations of current datasets 
point to physics beyond the minimal HZ model with stan- 
dard neutrino physics. In particular, allowing for either a 
scalar spectral index different from unity or an additional 
number of relativistic degrees of freedom, produces a sig- 
nificant increase in the goodness of fit with respect to the 
minimal HZ model. Information beyond that depends on 
the dataset used. 

For a dataset with ACT (and not SPT), there is no 
significant increase in the goodness of fit by increasing 
the complexity of the perturbation spectra by allowing 
for running of the scalar perturbations and/or a tensor 
component. Concerning neutrino physics, there is no sig- 
nificant further increase in the goodness of fit by allowing 
for a nonzero neutrino mass. For the dataset with ACT, 
we can only say that a model with a scalar spectral in- 
dex different than unity is a much better fit than the HZ 
model, and a model with a non-standard number of neu- 
trinos is also a much better fit. In the sections below we 
will quantify these statements. 

The situation is much different if we examine a dataset 
with SPT (and not ACT). 3 There are a significant in- 
creases in goodness of fit allowing either a running of the 
scalar spectral index (and not much increase in good- 
ness of fit just by allowing a tensor component) on top 
of the HZ+n s model, or a nonzero neutrino mass to the 
HZ model. Again, in the sections below we will quantify 
these statements. 

Therefore, we conclude that a cosmological dataset in- 
cluding SPT suggests either a more complex perturbation 
spectrum than simply a scalar spectral index different 
than unity, or some other new physics such as a modi- 
fication of the number of relativistic degrees of freedom. 
The data do not prefer one approach over the other. 

The papers is organized as follows: In the next section 
we review the pertinent features of slow-roll inflation. In 
Sec. |III| we review how neutrinos (or other light species) 



affect the CMB anisotropics. In Sec. [TV] we discuss our 
data analysis method and the datasets examined. Sec- 
tion [V] presents our results for cosmological parameters 
and the maximum likelihood for various models. In Sec. 



1 We note that a scalar spectral index of unity is possible within 
slow-roll inflation [TflfTB] . 

2 Although we parameterize the additional relativistic degrees of 
freedom as a contribution to the effective number of light neu- 
trinos, N e ff, of course the new relativistic species need not be 
neutrinos. 



3 Ref. 27 concluded that datasets including either ACT or SPT 
give a consistent picture for cosmological parameters, as long as 
YLZ+ns and HZ+n s + iV e ff are concerned. We find however that 
the answer to the question of whether cosmological data points 
to physics beyond an HZ +ns model strongly depends on the 
choice of whether the dataset includes ACT or SPT. 
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Vl| we discuss implications for physics beyond the HZ 



model for neutrino physics and for inflation. Section VII 
contains our conclusions. 



II. SLOW-ROLL INFLATION AND CMB 
ANISOTROPIES 

As mentioned in the introduction, we will work un- 
der the hypothesis that the adiabatic perturbations orig- 
inated within the single-field, slow-roll framework of in- 
flation. It should be kept in mind that if future exper- 
iments do not detect isocurvature modes or large non- 
Gaussianity it will not be possible to distinguish directly 
the inflaton contribution from the, e.g., curvaton con- 
tribution, see Ref. [28]. On the other hand, a detection 
of a significant amount of tensor modes through CMB 
anisotropics will disfavor curvaton-like models as they 
tend to generate a negligible tensor contribution. 

Within the single-field slow-roll paradigm, many spe- 
cific models for inflation have been proposed. We limit 
ourselves here to models with "normal" gravity (i.e., gen- 
eral relativity) and a single order parameter for the vac- 
uum, described by a canonical scalar field </>, the inflaton, 
with Lagrangian 
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The equations of motion for the spacetime are given by 
the Friedmann Equations, which for a homogeneous field 
6 are 
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The equation of motion for the field <f> is 
<f> + 3H4> + V (4>) = 0. 
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We have assumed a flat Friedmann-Robertson- Walker 
metric g^ v = diag(l, — a 2 , — a 2 — a 2 ), where a(t) is the 
scale factor of the universe. Inflation is defined to be a 
period of accelerated expansion, a > 0. If the field evo- 
lution is monotonic in time, we can write the scale factor 
a (</>) and Hubble parameter H (</>) as functions of the 
field (j) rather than time, i.e., we define all of our physi- 
cal parameters along the trajectory in phase space ((f)) 
corresponding to the classical solution to the equations 
of motion. Equations Q and ([3| can then be re- written 
exactly in the Hamilton-Jacobi form 
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These are completely equivalent to the second-order 
equation of motion. The second of the above equations 



is referred to as the Hamilton-Jacobi equation, and can 
be written in the useful form 
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where e is defined to be 
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The physical meaning of e((j>) can be seen by expressing 
a m Eq. (U) as 
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so that the condition for inflation, (a /a) > 0, is equiva- 
lent to e < 1. The scale factor is given by 
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where the number of e-folds N is 
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Most simple inflation models satisfy the slow-roll ap- 
proximation, which is the assumption that the evolution 
of the field is dominated by the drag from the cosmolog- 
ical expansion, so that </> ~ and (j) ~ —V/3H. The 
equation of state of the scalar field is dominated by the 
potential, so that p ~ — p, and the expansion rate is 
approximately H 2 c± SttV '((j)) / / 3mp 1 . The slow roll ap- 
proximation is consistent if both the slope and curvature 
of the potential are small, V V" <C V (in units of the 
Planck mass rapi). In this case the parameter e can be 
expressed in terms of the potential as 
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We will also define a second "slow-roll parameter" r] by 
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Slow roll is then a consistent approximation for e, ?|<Cl. 

Perturbations created during inflation are of two types: 
scalar (or curvature) perturbations, which couple to 
the stress-energy of matter in the universe and form 
the "seeds" for structure formation, and tensor, or 
gravitational- wave perturbations, which do not couple 
to matter. Both scalar and tensor perturbations con- 
tribute to CMB anisotropics. Scalar fluctuations can also 
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be interpreted as fluctuations in the density of the mat- 
ter in the universe. Scalar fluctuations can be quantita- 
tively characterized by the comoving curvature perturba- 
tion Pjz. As long as slow roll is attained, the curvature 
(scalar) perturbation at horizon crossing can be shown 
to be [T] 



dV2 

n 



(k) 





' H 


1 




^ k=aH 


mp i * 


s / / 7re_ 


k=a,H 



(12) 



The fluctuation power spectrum is, in general, a function 
of wavenumber /c, and is evaluated when a given mode 
crosses outside the horizon during inflation, k = aH. 
Outside the horizon, modes do not evolve, so the ampli- 
tude of the mode when it crosses back inside the horizon 
during a later radiation- or matter-dominated epoch is 
just its value when it left the horizon during inflation. 
Instead of specifying the fluctuation amplitude directly 
as a function of /c, it is convenient to specify it as a func- 
tion of the number of e-folds N before the end of inflation 
at which a mode crossed outside the horizon. 

The scalar spectral index ns for Pn is defined by 
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so that a scale- invariant spectrum, in which modes have 
constant amplitude at horizon crossing, is characterized 
by n s 1. 

To lowest order in slow roll, the power spectrum of 
tensor fluctuation modes and the corresponding tensor 
spectral index is given by pQ 
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The ratio of tensor-to-scalar modes is then Pt/Pu = 16e, 
so that tensor modes are negligible for e < 1. In the 
limit of slow roll, the spectral indices ns and tit vary 
slowly or not at all with scale. We can write the spectral 
indices ns and nr to lowest order in terms of the slow-roll 
parameters e and n as 



n s ~ 1 - 4e + 2r?, 
riT — — 2e. 
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The tensor /scalar ratio is frequently expressed as a 
quantity r, which is conventionally normalized as 



16e 



Pt 



(16) 



The tensor spectral index is not an independent param- 
eter, but is proportional to the tensor/scalar ratio, given 
to lowest order in slow roll by - — 2e = — r/8. A 
given inflation model can therefore be described to low- 
est order in slow roll by three independent parameters: 
P n , Pt, and n S - 



Deviations from a simple power-law spectrum of per- 
turbations are higher order in the slow-roll parameters, 
and thus serve as a test of the consistency of the slow-roll 
approximation. Scale dependence in the observables cor- 
responds to scale dependence in the associated slow-roll 
parameter, and can be quantified in terms of the infinite 
hierarchy of inflationary flow equations [29] , 
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The higher-order flow parameters are defined by 
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where the prime denotes derivatives with respect to 
scalar field (j). It is then straightforward to calculate the 
scale-dependence of the spectral index by relating the 
wavenumber k to the number of e-folds AT, 
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Since the running depends on higher-order flow param- 
eters than the spectral index itself, it is an independent 
parameter, even in slow-roll inflation models. In typical 
single-field inflation models, the running of the spectral 
index is negligible, so a detection of scale dependence in 
the spectral index would rule out a large class of viable 
single-field inflation models, and would therefore be a 
powerful probe of inflationary physics. 



III. NEUTRINOS AND CMB ANISOTROPIES 

In what follows we we examine the possibility of new 
neutrino physics as an alternative to extending the com- 
plexity of primordial perturbations. 

One direction for new neutrino physics is a change in 
the effective number of relativistic degrees of freedom, 
Aeff , that defines the physical energy density in relativis- 
tic particles p ra d, defined by 
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where p 1 is the energy density of the CMB photons and 
N e ff is the effective number of light neutrino species. In 
the standard scenario, assuming three active massless 
neutrino species with standard electroweak interactions 
and the present CMB temperature of T 7 = 2.726K (see, 
e.g., Ref. [30 J, the expected value is N e ^ = 3.046. This is 
slightly larger than 3 because of non-instantaneous neu- 
trino decoupling (see, e.g., Ref. [31]). As mentioned pre- 
viously, any new species that is relativistic around recom- 
bination will contribute to iV e ff, whether it is a neutrino 
species or not. The exact contribution of a new relativis- 
tic species will depend on the number of spin degrees of 
freedom, whether the new species is a boson or fermion, 
and the temperature of decoupling of the new species. 

We also consider the possibility of a mass for one or 
more of the three known active neutrino species. The 
present contribution to the overall energy density is given 

by 



=1,2,3 



92.5 eV 
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where are the masses of the three neutrino mass eigen- 
states. 

A change in neutrino physics can have important impli- 
cations for interpretation of inflationary parameters from 
CMB anisotropics, see Refs. [32H34] . For example, vary- 
ing iVeff can have an impact on determination of ns and 
its running, since it changes both the position of the CMB 
peaks in the angular spectrum and the structure of the 
"damping tail" at very large multipoles (see Ref. [TT]). 
In general, a higher N e ^ can put higher values of ns in 
better agreement with the data, i.e., there is a positive 
correlation between the two parameters. 

Masses for neutrinos also have important implications 
for interpretation of inflationary parameters from CMB 
anisotropics. Massive neutrinos damp the dark- matter 
fluctuations on scales below the horizon when they be- 
come nonrelativistic (see e.g., [35]). Neutrinos with 
masses m v < 0.3 eV are relativistic at recombination 
and affect the CMB anisotropy mainly through gravita- 
tional lensing, while neutrinos with larger masses slightly 
increase the CMB small-scale anisotropy by damping the 
gravitational potential at recombination. The final result 
is a small anti-correlation with n#, i.e., larger neutrino 
masses shift the constraints on ns to smaller values. 



IV. DATA ANALYSIS METHOD 

The analysis method we adopt is based on the publicly 
available Monte Carlo Markov Chain (MCMC) package 
cosmomc [36] with a convergence diagnostic done through 
the Gelman and Rubin statistic. 

We sample the following four-dimensional standard set 
of cosmological parameters, adopting flat priors on them: 
the baryon and cold dark matter densities and Q c , the 
angular size of the sound horizon at decoupling 9, and the 
optical depth to reionization r. 



As discussed in a separate section, we will also vary 
the relativistic number of degrees of freedom parameter 
iV e ff and the total neutrino mass Sm^. The standard 
three- neutrino framework predicts N e s — 3.046, while 
oscillation neutrino experiments place a lower bound 
Sm, > 0.05 eV [37]. 

For the inflationary parameters we consider the scalar 
spectral index ns and its running n mn , the overall nor- 
malization of the spectrum As at k = 0.002 Mpc -1 and 
the amplitude of the tensor modes relative to the scalar, 
r = A T /A s , again at k = 0.002 Mpc -1 . 

We consider purely adiabatic initial conditions and we 
impose spatial flatness. 

We analyze the following set of CMB data: WMAP9 
[19], SPT [20], and ACT [21], including measurements up 
to a maximum multipole number of / max ~ 3000. For all 
these experiments we make use of the publicly available 
codes and data. For the ACT experiment we use the 
"lite" version of the likelihood [38] . 

We also consider the effect of including additional 
datasets to the basic datasets just described. Consis- 
tently with the measurements of HST [25 , we consider 
a Gaussian prior on the Hubble constant Ho = 73.8 ± 
2.4 km s -1 Mpc -1 . We also include information from 
measurements of baryonic acoustic oscillations (BAO) 
from galaxy surveys. Here, we follow the approach pre- 
sented in Ref. [19] combining three datasets: SDSS-DR7 
[22], SDSS-DR9 [23] and WiggleZ [24]. 

Since, as we see in the next section, the ACT and the 
SPT datasets are providing significantly different conclu- 
sions on inflationary parameters, we will include them 
separately. In what follows we will consider two combi- 
nations of datasets. We refer to an analysis using the 
WMAP9 + ACT + HST + BAO datasets as the "ex- 
tended ACT" (eACT) dataset and to an analysis with 
the WMAP9 + SPT + HST + BAO datasets as the "ex- 
tended SPT" (eSPT) dataset. 

We use the Markov chains obtained from CosmoMC 
to reconstruct the posterior distributions of each of the 
model parameters. In the tables, we present our results in 
the form of the 68% credible interval for each parameter, 
i.e., the interval symmetric around the mean containing 
68% of the total posterior probability. We make an ex- 
ception to this rule in those cases where the posterior 
probability is not vanishingly small at the edge of the 
prior range; this happens in particular around r = and 
^2 m v = 0. In this case we adopt the following rule: if 
the maximum of the posterior distribution is clearly dis- 
tinguished from zero, we quote the 68% interval as above; 
otherwise, we quote a 95% upper limit. 

We also use our Markov chains to recover the maxi- 
mum likelihood (i.e., minimum y 2 ) parameter values. We 
use the minimum x 2 values estimated from the chains to 
perform an approximate model comparison by comput- 
ing the likelihood ratio (actually, equivalently, the differ- 
ence in x 2 ) between models. As a rule of thumb, given 
two models M.\ and M2, where the latter reduces to 
the former for a particular choice of parameter values (in 
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which case the two models are said to be "nested"), we 
say that the data show preference for M2 over Mi when 
the absolute value of Ax 2 = XminO^) - Xmm( M i) is 
larger than the number of additional parameters in the 
extended model. 

We note however that MCMC methods are usually op- 
timized to sample the full posterior distribution around 
the region of maximum probablity, and not to recover the 
exact value and position of the maximum likelihood. The 
precision to be associated with our estimate of the mini- 
mum x 2 can De evaluated by computing the probability 
of finding in the chains a sample having a x 2 within S(x 2 ) 
from the actual minimum (see e.g. [39]). This probabil- 
ity depends on the dimensionality of the parameter space 
and on the number of independent samples in the chains. 
We let our chains run until we can claim a 95% probabil- 
ity of having found the best-fit model with an uncertainty 
% 2 ) < I- 



V. EXTENSIONS OF THE HZ MODEL AND 
EACT AND ESPT 

As a first step in our analysis we evaluate the compat- 
ibility of current cosmological datasets eSPT and eACT 
with a simple reference model, which we choose to be 
the Harrison-Zel'dovich (HZ) model with ns = 1, r = 0, 
n run = 0, m v = 0, and 7V" e ff = 3.046. We then consider 
extensions of this model involving more complex pertur- 
bation spectra, with various combinations of ns 7^ 1, 
r > 0, and n run =^ 0. Then we examine extensions of 
the HZ model with nonstandard neutrino physics with 
combinations of m v ^ and/or 7V eff ^ 3.046. 



A. Extensions of the perturbation sector 

The results of our analysis with regard to perturbation 
spectra is reported in Table [IJ As stated in the previous 
section, we analyze the eACT and eSPT datasets and we 
consider different cases for primordial perturbations and 
compare them with the reference HZ model. In all models 
analyzed in this section we assume massless neutrinos 
and 7V eff = 3.046. 

As we can see from the table, both for the eACT and 
eSPT datasets, models with ns ^ 1 are highly favored 
over the HZ reference model. 

For the eSPT dataset, allowing one additional param- 
eter, ns, to vary results in change in x 2 °f ^X 2 = 
(— 21og£) — (— 2\og Chz) = —28.7. The one-dimensional 
probability distribution for ns with the eSPT dataset 
is shown in Fig. [I] For the eACT dataset, allowing 
one additional parameter, ns, to vary results in Ax 2 = 
(— 21og£) — (— 2\og Chz) = —9.7. The one-dimensional 
probability distribution for ns with the eACT dataset is 
also shown in Fig. [I] 

If we allow other parameters describing the perturba- 
tion spectra to vary, such as n mn and r, there are differ- 



ent indications from the different datasets. Let us first 
consider the eSPT dataset. 

The natural parameter space for constraining sim- 
ple slow-roll inflation models is to include the the ten- 
sor/scalar ratio r in addition to spectral tilt ns- Two- 
dimensional contours for ns vs. r are shown in Fig. [2j 
along with the predictions of three simple slow-roll mod- 
els. For the eSPT dataset, allowing r to vary in addi- 
tion to allowing ns to vary results in a very marginal 
decrease in x 2 of —0.9 compared to a model just allow- 
ing ns to vary. Hence, the data do not seem to call 
for the additional variable r. However, the situation is 
quite different if we allow a running of the scalar spec- 
tral index, n run 7^ 0, either keeping r = or allowing r 
to vary. Adding one additional parameter, n run , results 
in Ax 2 = —36.6 compared to the reference HZ model, 
which corresponds to a change in x 2 of —7.9 compared 
to the HZ +ns model. If we allow both r and n run to 
vary (in addition to allowing ns to vary) there is a gain 
of Ax 2 = —37.5 compared to the reference HZ model, or 
a change in x 2 of —8.8 compared to the HZ +ns model. 
The eSPT dataset strongly prefers a running of the scalar 
spectral index. The one-dimensional probability distri- 
butions for ns, r and n run with the eSPT dataset are 
shown in Fig. [I] Two-dimensional contours of r vs. ns, 
n run vs. ns and n run vs. r are also shown in Fig. [I] 

The eACT dataset also prefers a scalar spectral index 
different from unity. Recall that adding one additional 
parameter ns results in a decrease in x 2 compared to the 
reference HZ model of Ax 2 = — 9.7. If we then allow one 
additional parameter, either r or n mn , there is only a very 
marginal change in x 2 beyond the HZ + ns model. Even 
allowing both additional parameters n mn and r again re- 
sults in a very marginal decrease in x 2 at the expense of 
two additional parameters. The one-dimensional prob- 
ability distributions for ns, r and n run with the eACT 
dataset are also shown in Fig.[l] and the two-dimensional 
contours of r vs. ns, n run vs. ns and n run vs. r are also 
shown in Fig. [T] 

We summarize our findings with respect to n s and n mn 
in Fig. [3J where we compare the constraints on these 
parameters for the different model/dataset combinations 
considered in the paper. It is clear from this figure that 
the tension between the two datasets increases when the 
model complexity is also increased. Moreover, as dis- 
cussed above in the context of the goodness-of-fit of the 
various models, we also notice that the results of param- 
eter estimation from eACT are more stable, with respect 
to eSPT, to the increase of the complexity of the model. 

Our conclusion is that the eSPT and eACT datasets 
are not consistent, as long as inflation- motivated exten- 
sions to the minimal model are concerned. While both 
call for a scalar spectral index different than unity, the 
eSPT dataset seems to be better described by a more 
complicated perturbation spectrum than just a scalar 
spectrum of constant spectral index. On the other hand, 
the eACT dataset seems to be well described by a con- 
stant scalar spectral index (slightly less than unity) , and 
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TABLE I: Augmenting the minimal Harrison- Zel'dovich cosmological model through inflationary parameters. Listed are pos- 
terior means for the cosmological parameters from the indicated datasets (errors refer to 68% confidence intervals, unless 
otherwise stated). 







Reference Model 




Inflation-Motivated Extensions 




Dataset 


Parameter 
















HZ 


HZ ±n s 


HZ ±n s ± r 


HZ +n s ± n run 


HZ +ns ± r ± n run 




100 Q b h 2 


2.331 ± 0.025 


2.225 ±0.032 


2.228 ±0.032 


2.236 ± 0.031 


2.272 ±0.036 




n c h 2 


0.1148 ±0.0017 


0.1167 ±0.0018 


0.1166 ±0.0018 


0.1180 ±0.0019 


0.1178 ±0.0018 




100 


1.0430 ±0.0009 


1.0419 ±0.0009 


1.0419 ±0.0010 


1.0422 ±0.0009 


1.0424 ±0.0009 




log[10 10 A s ] 


3.12 ±0.03 


3.21 ±0.03 


3.20 ±0.03 


3.14 ±0.04 


3.04 ±0.07 


eSPT 


T 


0.096 ±0.013 


0.078 ±0.012 


0.077 ±0.012 


0.090 ±0.014 


0.095 ±0.015 


n s 


= 1 


0.959 ±0.008 


0.962 ± 0.008 


1.037 ±0.029 


1.107 ±0.045 




r 








< 0.12 (d) 





0.28 ±0.16 




?T/run 


= 


= 


= 


-0.029 ±0.011 


-0.051 ±0.015 




Ho (a) 


71.33 ±0.65 


69.33 ± 0.74 


69.42 ± 0.76 


69.08 ±0.76 


69.51 ±0.78 




-21og£ (b) 


7653.4 


7624.7 


7625.6 


7616.8 


7615.9 




Ax 2 (c) 





-28.7 


-27.8 


-36.6 


-37.5 




100 Q b h 2 


2.356 ± 0.027 


2.282 ±0.035 


2.290 ±0.037 


2.283 ± 0.035 


2.302 ± 0.038 




n c h 2 


0.1163 ±0.0021 


0.1165 ±0.0021 


0.1162 ±0.0021 


0.1166 ±0.0021 


0.1167 ±0.0022 




100 


1.0416 ±0.0016 


1.0399 ±0.0018 


1.0399 ±0.0017 


1.0400 ±0.0017 


1.0403 ±0.0018 




log[10 10 A s ] 


3.14 ±0.03 


3.19 ±0.03 


3.18 ±0.03 


3.19 ±0.04 


3.13 ±0.05 


eAC 1 


r 


0.102 ±0.014 


0.090 ±0.014 


0.089 ±0.013 


0.092 ±0.015 


0.094 ±0.015 


n s 


= 1 


0.971 ±0.009 


0.976 ± 0.009 


0.978 ± 0.031 


1.016 ±0.042 




r 








< 0.18 (d) 





< 0.34 (d) 




Tlrun 


= 


= 


= 


-0.003 ±0.011 


-0.014 ±0.014 




H <«> 


70.50 ±0.71 


69.24 ±0.83 


69.43 ± 0.83 


69.24 ±0.81 


69.47 ±0.83 




-21og£ (6) 


7617.9 


7608.2 


7608.4 


7608.3 


7608.7 




Ax 2 (c) 


= 


-9.7 


-9.5 


-9.6 


-9.2 



a km s -1 Mpc -1 

fe When comparing to the x 2 values reported e.g. in the WMAP9 
paper 19 , it should be taken into account that we use a pixel based 
likelihood at low Is instead than the Gibbs-based likelihood. 

C A X 2 = (-2 log/:) - (-2 log £ H z) 

d 95% c.l. 



does not seem to require additional complexity. 



B. Extensions of the neutrino sector 

We now repeat the analysis presented in previous sub- 
section but now considering the possibility of an extra 
effective neutrino number and including neutrino masses. 

The constraints from the eSPT and eACT dataset are 
in Table [TTJ For both datasets, adding the additional pa- 
rameter N e f[ greatly improves the fit. In fact, allowing 
N e ff improves the fits of both eSPT and eACT by about 
as much as allowing the spectral index to vary from unity. 



However, allowing the neutrino mass to vary, we again 
again obtain different indications from the two datasets. 
For the SPT dataset, adding a neutrino mass improves 
the x 2 by —8.1 if N e ^ is kept fixed and by —8.8 if it is 
allowed to vary. For the ACT dataset, on the contrary, 
the goodness of fit improves only marginally (at the price 
of one additional parameter) by allowing a non-zero neu- 
trino mass, independently of whether N e s is fixed or not. 

In Fig. [5] we compare the constraints on N e ^ and m v 
for the different model/dataset combinations considered 
in the paper. Again we see the same trend observed 
in the case of the spectrum parameters, namely that the 
values estimated from the two datasets tend to diverge as 
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FIG. 1: One- and two-dimensional posterior probabilities for ns, r, and n vun - Upper panel: One-dimensional parameter 
posteriors for the models considered in the text, using the eSPT (top row) and eACT (bottom row) datasest. Lower panel: 
Two-dimensional posteriors for the HZ+n s +r+running case. Dark- and light-shaded regions correspond to 68 and 95% credible 
intervals, respect ivey. 



new parameters are added, and that the values estimated 
from eACT are more stable than those estimated from 
eSPT when the complexity of the model is increased. 



VI. DIRECTIONS FOR NEW PHYSICS 

Using the two data combinations described here, mod- 
els with either primordial perturbations beyond the HZ 
model or additional light degrees of freedom provide a 
much better fit than the HZ model. Thus, cosmological 
data point to some interesting new physics. Unfortu- 
nately, the direction is unclear. 



Probably the most dramatic explanation would be ad- 
ditional light degrees of freedom: It would be very sur- 
prising if there is a new light species beyond the standard 
model of particle physics (as we have emphasized, it need 
not be extra neutrino species, although we parameterize 
them as such). 

For the eACT dataset, just adding a tilt to the scalar 
spectrum seems to be all that is demanded of the data. 
This would tell us something about inflation, but there 
are a large number of inflation models that can give a 
slightly red spectrum. 

For the eSPT dataset however, the data seems to de- 
mand more than simply a tilt to the scalar spectrum. A 
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FIG. 3: Comparing the constraints on n s (top panel) and 
n run (bottom panel) for different model/datasets combina- 
tions. The solid and dashed bars denote 1- and 2-a con- 
straints, respectively . 



much improved fit can be obtained by allowing the possi- 
bility of a large running of the scalar spectrum. The run- 
ning could be so large as to have a large impact in infla- 
tion model building and call in doubt the simple slow-roll 
approximation. Alternatively, as data seem to indicate 
a non-power-law scalar spectral index with a very large 
variation of the spectral index, one might invoke models 
where the flattening of the inflaton potential is obtained 
through the inclusion of large quantum corrections in the 
mass parameter [40] which result in large variation of the 



spectral index with the scale. Another class of models 
which allow for a large negative running are models in 
which inflation occurs near an inflection point of the po- 
tential, where the third derivative V'" of the potential 
is substantial, and the higher-order slow roll parameter 
2 A is comparable to the lower-order parameters e and n. 
Inflection point inflation models have been argued, e.g. 
in Ref. [41] , to be characteristic of inflation on the string 
landscape. 



VII. CONCLUSIONS 

We analyzed the recently released Atacama Cosmol- 
ogy Telescope (ACT) and South Pole Telescope (SPT) 
data in combination with the Wilkinson Microwave 
Anisotropy Probe 9-year data (WMAP9), the Sloan Dig- 
ital Sky Survey Data Release 9, the WiggleZ large-scale 
structure data, and the Hubble Space Telescope deter- 
mination of the Hubble parameter (HST). We tested 
these data against two cosmological scenarios: (1) a 
scale-invariant, purely scalar "Harrison-Zel'dovich" (HZ) 
power spectrum with the addition of parameters moti- 
vated by inflationary cosmology, tilt ns, nonzero ten- 
sor/scalar ratio r, and running of the spectral index ns, 
and (2) the HZ power spectrum with a nonstandard ef- 
fective neutrino number N e ff and/or neutrino mass m v . 
We find that both the extended ACT data (eACT) and 
the extended SPT data (eSPT) favor extensions to the 
simple HZ model to at least 95% confidence. 

In the case of the inflation- motivated extensions to HZ, 
both eACT and eSPT favor a deviation from a scale- 
invariant power spectrum with "red" tilt, ns < 1, and 
neither show any evidence for a nonzero tensor/scalar 
ratio. The eACT data are consistent with negligible 
running of the spectral index, as predicted by simple 
slow-roll inflation models. The eACT data are consis- 
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TABLE II: Augmenting the minimal Harrison-Zel'dovich cosmological model through new neutrino physics. Listed are posterior 
means for the cosmological parameters from the indicated datasets (errors refer to 68% credible intervals, unless otherwise 
stated) . 



Dataset 


Parameter 


Reference Model 


Neutrino-Motivated Extensions 






HZ 


TT V _|_ A! 




ti y _|_ at- _i_ 




100ftb/i 2 


2.331 ±0.025 


o o i i _l n no a 
Z.oll ± U.Uz4 


o oon _l n no a 
z.ooU ± U.Uz4 


o ooo _i_ n noT 

z.332 ± U.U37 




Q c h 2 


0.1148 ±0.0017 


n i on/1 _l n nncv 

U.1394 ± U.UU57 


n 1 1 nn _l n nnoo 

U.11UU ± U.UUzo 


n 1 oi k _i_ n nncv 

U.lolo ± U.UU57 




100 


1.0430 ±0.0009 


i r\Ar\A _i_ n nmn 
1.U4U4 ± U.UU1U 


i.u4o4 ± u.uuuy 


1 n /i 1 o _l n nm 1 
1.U41Z ± U.UU11 




\og[10 10 A s ] 


3.12 ±0.03 


Q 1 K _i_ n HQ 
o.lo zb U.Uo 


o.lz zt U.Uo 


q i/| i n no 
o.l4 zh U.Uo 


eSPT 


T 


0.096 ±0.013 


n no c _i_ n m o 
U.UoO ± U.Ulz 


n i no _i_ n ni /i 
U.lUo zt U.U14 


n nnpc _l n m /i 
U.Uyo zt U.U14 


iVeff 


= 3.046 


4. ZD ± U.zO 


— O.U4D 


A A K -\- fl Q9 
4.40 zt U.oz 









= u 


U.oy ± U.14 


n n« _i_ n co 
u.yo zt U.oo 




Ho 


71.33 ±0.65 


75.5 ±1.1 


69.82 ±0.76 


74.0 ± 1.2 




-21og£ (c) 


7653.4 


/ ozo.y 


IRA K Q 


/ Ol / .1 




A X 2 (d) 





— Z / .0 


— O.l 


— OD.O 




100Qb/i 2 


2.356 ± 0.027 


o ooo j_ n non 
z.ooz ± U.Uzy 


o oro _i_ n non 
z.ooo ± U.Uzy 


o qq7 j_ n non 
z.oo / ± U.Uzy 




Q c h 2 


0.1163 ±0.0021 


n 1 oi o _i_ n nn ^7 

0.1318 ± 0.0057 


n 1 1 cr/^ _i_ n nnoi 

U.lloo ± U.UUzl 


n 1 on/ 1 _i_ n nntr^ 

0.1z9o ± 0.U057 




100 6 


1.0416 ±0.0016 


i nqoo -4- D 0090 
j-.uooz ztz u.uuzu 


i 0491 4- 001 

1 .U^Z ± ZC U.UU1U 


i noo7 _i_ n 0090 
i.uooi zn u.uuzu 




\og[10 10 A s ] 


3.14 ±0.03 


3.16 ±0.03 


3.13 ±0.03 


3.15 ±0.03 


eACT 


T 


0.102 ±0.014 


0.097 ±0.014 


0.105 ±0.015 


0.099 ±0.014 




iVeff 


= 3.046 


3.88 ±0.28 


= 3.046 


3.80 ±0.28 









= 


0.24 ±0.15 


< 0.46 (e) 




H <*> 


70.50 ±0.71 


73.2 ±1.1 


69.82 ±0.79 


72.4 ± 1.2 




-2 log/: (c) 


7617.9 


7609.7 


7616.7 


7609.2 




A X 2 (d) 





-8.2 


-1.2 


-8.7 



a eV 

5 km s -1 Mpc -1 

c When comparing to the % 2 values reported e.g. in the WMAP9 
paper 19 , it should be taken into account that we use a pixel based 
likelihood at low Is instead than the gibbs-based likelihood. 

d A x 2 = (-2 log/:) - (-21og£^ z ) 

e 95% c.l. 



tent at the 95% confidence level with simple chaotic 
inflation V (0) = m 2 (/> 2 , and with power-law inflation, 
V (0) oc exp(aV/i), as well as "small-field" models pre- 
dicting negligible tensors and ns < 1. The eSPT data, 
however, are inconsistent with a purely power-law power 
spectrum, favoring negative running of the spectral index 
n ron = —0.029 ± 0.011 in the case with a prior of r = 0, 
and n Tun = —0.051 ± 0.015 in the case where r ^ is 
allowed. While the eSPT data are not in disagreement 
with the most general possible single-field inflation mod- 
els, they are in significant conflict with slow-roll models 
predicting n run <C n^. The eACT data are consistent 
with such models. 

In the case of extensions to HZ involving additional 
light degrees of freedom, eACT and eSPT again produce 



qualitatively different constraints. Both the eACT and 
eSPT data favor additional light degrees of freedom, with 
N eff = 3.88 ±0.28 for eACT, and N eff = 4.26 ±0.26 for 
eSPT (with a prior of m v = 0). The eACT and eSPT 
data differ, however, with respect to nonzero neutrino 
masses. The eACT data are consistent at 95% with zero 
neutrino mass, with 52 m v = 0.24±0.15 eV (with a prior 
of N ef f = 3.04), and ^m, < 0.46 eV (with N eff ^ 
3.04). The eSPT data favor nonzero neutrino mass, with 
52 m u = 0.39 ± 0.14 eV (with a prior of N eff = 3.04), 
and 52m u = 0.96 ± 0.53 eV (with N eff ^ 3.04). 

In either scenario, HZ ± inflation or HZ ± neutrinos, 
considering the ACT and SPT data separately results 
in qualitatively different conclusions about extensions to 
a standard scale- invariant A±Cold Dark Matter concor- 
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dance cosmology, a tension which is not evident when 
considering combined constraints from ACT and SPT. 
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